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We investigate spherically symmetric, static spacetimes in Eddington-inspired Born-Infeld gravity 
coupled to Born-Infeld electrodynamics. The two constants, 6^ and k which parametrise the Born- 
Infeld structures in the electrodynamics (matter) and gravity sectors, characterise the features of our 
analytical solutions. Black holes or naked singularities are found to arise, depending on the values 
of 6^ and k, as well as charge and mass. Several such solutions are classified and understood through 
the analysis of the associated metric functions for fixed k, varying and vice-versa. Further, we 
also compare the new metric functions with those for the known —>■ oo (Maxwell) and the re —>■ 0 
(geonic black hole) cases. Interestingly, for a particular relation between these two parameters, 
b^ = l/4re, re > 0, we obtain a solution resembling the well-known Reissner-Nordstrom line element, 
albeit some modifications. Using this particular solution as the background spacetime, we study 
null geodesics for Born-Infeld photons and also, gravitational leasing. Among interesting features 
we note (i) an increase in the radius of the photon sphere with increasing re and (ii) a net positive 
contribution in the leading order correction term involving re, in the weak leasing formula for the 
deflection angle. We also investigate the effective potential and light propagation for various other 
solutions through numerics and plots. In summary, our work is the first attempt towards figuring 
out how Born-Infeld structures in both the matter and gravity sectors can influence the nature and 
character of resulting gravitational fields. 
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I. INTRODUCTION 


General Relativity (GR) is largely successful as a classical theory of gravity. It has been tested through experiments 
and observations to very high precision, in vacuum. However, it is not entirely certain that GR is the correct theory 
of gravity inside a matter distribution or, in the strong field regime. Moreover, it possesses many unsolved puzzles, 
both theoretical and observational. This has prompted many researchers to actively pursue various modified theories 
of gravity. For example, one motivation for constructing modified theories of gravity is to look for the possibility of 
avoiding the singularity problem in GR, eg. the occurrence of a big bang in cosmology or black holes in an astrophysical 
context. In a classical metric theory of gravity, under very general and physical assumptions, these singularities are 
inevitable, as proved through the Hawking-Penrose theorems [l|. However, it is possible to obtain non-singular 
geometries as solutions in an alternative theory. For example, non-singular Friedman-Robertson-Walker (FRW) type 
cosmologies do exist in modified gravity theories. One can also argue that the singularity in the classical theory is 
expected to be resolved in quantum gravity. Unfortunately, though specific examples exist in the context of loop 
quantum cosmology, we are still far away from any general statements which may be possible when we have a fully 
consistent quantum description of gravity. 

In this article, we investigate one such alternate theory which is inspired from Eddington’s alternative formulation of 
GR. In this formulation, the y/—gR (Ricci scalar) in the Einstein-Hilbert action is replaced by y^—det{Rij) and the 
connection is chosen as the basic variable, instead of the metric. Thus, Eddington’s formulation is essentially an affine 
formulation Q. In vacuum, it is equivalent to GR and it describes de Sitter gravity. However, coupling of matter was 
never addressed in this formulation until very recen^. 

We are also aware of Born-Infeld electrodynamics Q wherein the divergence in the electric field at the location of 
the charge/current is removed by using an action different from that of Maxwell. The new field theory due to Born 
and Infeld includes the determinantal structures in the action, inspired by Eddington’s work. The virtue of this 
nonlinear electrodynamics in removing the singularity of the electric field, inspired Deser and Gibbons to construct a 
Born-Infeld gravity theory in the metric formulation Q. Later, Vollick Q introduced the Palatini formulation in Born- 
Infeld gravity and worked on various related aspects. He also introduced a non-trivial and somewhat artificial way 
of coupling matter in such a theory [^, Q . More recently, Banados and Ferreira Q have come up with a formulation 
where the matter coupling is different and simpler compared to Vollick’s proposal. We will focus here on the theory 
proposed in Q and refer to it as Eddington-inspired Born-Infeld (EiBI) gravity, for obvious reasons. Note that the 
EiBI theory has the feature that it reduces to GR, in vacuum. 

Interestingly, EiBI theory also falls within the class of bimetric theories of gravity (also called bi-gravity). The current 
bimetric theories have their origin in the seminal work of Isham, Salam and Strathdee Several articles have 
appeared in the last few years on various aspects of such bi-gravity theories. In 0, the authors have pointed out 
that the EiBI field equations can also be derived from an equivalent bi-gravity action. This action is closely related to 
a recently discovered family of unitary massive gra vity theories which are built as bi-gravity theories. Several others 
have contributed in this direction, in various ways [111 . [l2| . 

Let us now briefly recall Eddington-inspired Born-Infeld gravity. The central feature here is the existence of a physical 
metric which couples to matter and another auxiliary metric which is not used for matter couplings. One needs to 
solve for both metrics through the field equations. The action for the theory developed in Q, is given as: 


'S'B/(ff,r,'I') 


1 

Sttk 


d^3 


\J + kR^,j{T)\ Ai/ 'g 


Siviig, 'h) 


( 1 ) 


where A = We assume G = 1 and c = 1 throughout in the paper. Variation w.r.t F, done using the auxiliary 
metric 9^1, = -I- kR^^( 9) gives 


Qfiu — ~b eiR^i/(^g') 


( 2 ) 


Variation w.r.t g^v gives 


= Xy/^g^" - ( 3 ) 

In order to obtain solutions, we need to assume a g^i, and a with unknown functions, as well as a matter stress 
energy Thereafter, we write down the field equations and obtain solutions using some additional assumptions 

about the metric functions and the stress energy. 

Quite some work on various fronts have been carried out on various aspects of this theory, in the last few years. 
Astrophysical scenarios have been discussed in [l 3 l - l 2 ll| while cosmology has been dealt with i n fl^ [ 23 - 1 ^ . Among 
other topics, a domain wall brane in a higher dimensional generalisation, has been analysed in [ 30 (|. Generic features 
of the paradigm of matter-gravity couplings have been analysed in [ 3 l|. In [s^, authors have shown that EiBI 
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theory admits some nongravitating matter distribution, which is not allowed in GR. Some interesting cosmological 
and circularly symmetric solutions in 2 + 1 dimensions are shown in (3^ . Constraints on the EiBI parameter k 
are obtained from studies on compact stars El [13, tests in the solar system [^, astrophysical and cosmological 
observations [l^, and nuclear physics studies [Sg- In [s^, a major problem related to surface singularities has been 
noticed in the context of stellar physics. However, gravitational backreaction has been suggested as a cure to this 
problem in In [s^, the authors propose a modification to EiBI theory by taking its functional extension in 

a way similar to f{R) theory. Recently, the authors in use a different way of matter coupling, by taking the 
Kaluza ansatz for five dimensional EiBI action in a purely metric formulation and then compactify it using Kaluza’s 
procedure to get a four dimensional gravity coupled in a nonlinear way to the electromagnetic theory. 

Not much has however been done on black hole physics, or, broadly on the spherically symmetric, static solutions 
in this theory. It may be noted that the vacuum, spherically symmetric static solution is trivially same as the 
Schwarzschild de Sitter black hole. But, the electrovacuum solutions are expected to deviate from the usual Reissner- 
Nordstrom solution in GR. This has been shown in [^, [^ where the authors consider EiBI gravity coupled to a 

Maxwell electric field of a localized charge. They obtain the resulting spacetime geometries, and study its properties. 
The basic features of such spacetimes includes a singularity at the location of the charge which may or may not be 
covered by an event horizon. The strength of the electric field remains nonsingular as in Born-Infeld electrodynamics. 
However, this may not be the only solution because, in EiBI gravity, the matter coupling is nonlinear. The authors 
of [d^l have shown, in a different framework, that the central singularity could be replaced by a wormhole supported 
by the electric field. In [d^ , the author obtained a class of Lorentzian regular wormhole spacetimes supported by the 
quintessential matter which does not violate the weak or null energy cindition in EiBI gravity. 

In our present work, we couple EiBI gravity with Born-Infeld electrodynamics. Thus, we have Born-Infeld structures 
in both the gravity and matter sectors. We obtain some new classes of spherically symmetric static spacetimes and 
study their properties. The solutions we present include black holes and naked singularities. Earlier, a lot of work has 
indeed been done by considering nonlinear electrodynamics coupled to GR |ddl - l5l| . Some of them are motivated by 
string theory since Born-Infeld structures naturally arise in the low energy limit of open string theory Is^ . Our 
article is sectioned as follows. In Section H, we derive the form of the stress-energy tensor assuming an ansatz for the 
physical metric {gfiu)- We obtain the general expressions for the metric functions by solving the EiBI field equations 
with an ansatz for the auxiliary metric {q^u)- In Section HI, we list the solutions by scanning across regions in the 
parameter space (6^, k), where 6^ and k are the parameters in the electromagnetic and gravity sectors, respectively. 
In Section IV, we analyze, in detail, a particular class of solutions of the Reissner-Nordstrom-type for which 6^ = ^ 
(k > 0). Here we also discuss null geodesics, photon propagation, and gravitational lensing. Finally, in Section V, we 
investigate the effective potential and gravitational lensing for the ^ solutions. In Section VI, we summarize 
our results and indicate possibilities for future work. 


II. SPHERICALLY SYMMETRIC STATIC SPACETIMES DUE TO A CHARGED MASS 


In curved spacetime, the Lagrangian density for the Born-Infeld electromagnetic field theory is given by Q, 


C = 


47r 


1 - 


62 64 


( 4 ) 


where, F = and G = are two scalar quantities constructed from the components of the electro¬ 

magnetic field tensor {F^^y) and the dual field tensor {G^.v)- Here, 6 sets an upper limit on the electromagnetic field 
and, when 6 —>■ oo. Maxwell’s theory is recovered. The resulting energy-momentum tensor has the following general 
expression; 


rj. ^ _2 dC 


Att 


In our work here we restrict ourselves to an electrostatic scenario (i.e. = {(/), 0,0, 0}, 1/ = 0), in the background of 

a spherically symmetric, static spacetime. We assume the line element to be of the form, 

ds2 ^ + U{r)e^’^^^'>df^ + V{f)f^ {d0^ + sin^ (6) 


where, r is the radial coordinate, but the physical radial distance or, the area radius of each 2-sphere is r = r\JV (r). 
With this assumption, solving the equations of motion for the electromagnetic field, we obtain 

^ ( 7 ) 


d(j) 


dr 
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were, q is the electric charge. In the Minkowski limit of Eq. (jG]), the Eq. ([7]) resembles the known Born-Infeld electric 
field due to a point charge. The non-zero components of the energy-momentum tensor are given as: 


Ttt = 


b'^U ( 


47r 


Ef2 


- 1 e 


2V> 




47r 


Tee = 


62 


1 - 


Er2 


47r I y/y2^4 _|_ j\p. 


62 

and, Trhrh = — I 1 - 


Ef2 


47r I ^y2^4 _|_ q2lJj2 


Ef2, 


Vr^ sin2 6. 


( 8 ) 


The energy-momentum tensor is not traceless ( i.e. ^ 0 for any q/b ^ 0). The auxiliary line element is assumed 

to be of the form 


dsl = + e^''df^ + {dd^ + sin2 ddcj)^) (9) 

Here, the area radius of the auxiliary spacetime is f. It is worth noting that our choice of the physical and auxiliary 
line elements are opposite to the standard choices in l40l - l^ . While making this choice, we have exploited two 
facts: {i) the absence of differential expressions in the field equations in ((S]) which opens up the scope of algebraic 
manipulations, and (ii) the specific forms of the tt- and ff- components of the stress-energy tensor. This choice 
simplifies the situation to quite an extent. Using Eqs. dH),®, and ® in the field equations obtained from the ‘g’- 
variation [Eq. ®], and after some algebra , we obtain the expressions for U and V. We get a quadratic equation for 
V as. 


(1 - 4k62) - 2 (1 - 2Kb^) U -f 1 - ^ = 0 (10) 

There are two possible solutions for V from Eq. (ITUl) and consequently, two solutions for U. We choose the solution for 
U and V such that U —1, U —1 for large f and the physical spacetime is that found in the GR regime. Therefore, 
we have. 


V(f) 

U(f) 


1 — 2k62 



1 — 4k62 




1 -I- -2— 

r -I- {,2p4 


_ 4^ _ 2k62 


(1 — AKb‘^)\J 


1 -I- -2— 


4:Kq^ 


From the field equations obtained via T’-variation [Eq. ®], we get, 

Ke-^'>’Rtt{q)-U = -l, 


Ke 


-2u 


Rrr {q) 


and 


yf2 


-d?<" 


u = 1 
) + i' 


= r2, 


( 11 ) 

( 12 ) 


(13) 

(14) 

(15) 


where, Rttiq) = e 2i'-i-2V' ^ 2 |— ^,2 _|_ Rff{q) = ^ + — ip”. Here, primes denote derivatives 

with respect to r. From Eq. (I13p and Eq. da, we get Ip' + i>' = 0. So, without loss of any generality, we can assume 
V = —Ip. Thus, from the Eq. (USD , we arrive at, 


1 

g2v.(r) = 1 _ !_ + _L 
3k Kr 


V(r)r'^dr -|- 


C 


(16) 


where, C is a constant of integration. For, q = 0 {i.e. vacuum), U = V = 1 [Eqs. (flTll . dllD ] for any k, and the solution 
should be the Schwarzschild spacetime. Therefore, C is identified as related to the total mass, i.e., C = —2M. 
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III. CLASSIFICATION OF SPACETIMES 


A. General Solutions 


In this section, we obtain different classes of solutions by scanning the b'^ — k parameter space. For this purpose, 
we introduce the definition b"^ = ^- We note that, if k > 0, then a > 0, and if k < 0 then a < 0. The characteris¬ 
tics/forms of the solutions change over different ranges of a. These are-(i) a —> oo, {ii) oo > a > 1, (iii) a —>• 1, {iv) 
—oo < a < 1 (but 7 ^ 0), and (u) a —>■ —oo. We discuss the a = 1 case elaborately, as a special case, in a separate 
section. The remaining classes of solutions are analysed here. Firstly, we note that a —>■ ±oo cases are actually those 
belonging to the Maxwell limit, i.e., a coupling of EiBI gravity with a Maxwellian electric field. This has already 
been studied earlier (dol - ld^ . One can verify that the solution for this case (particularly for a —>■ oo), after rewriting 
in the Schwarzschild gauge, takes exactly the same form as that of [i^. Next, we obtain the solutions for the cases 
(ii) and (iv) below: 


First, we rewrite U, V as, 


U(r) = 
V(f) = 


2-a 


a 


2(1 — a) 2(1 — q) L 4K.q^il-a) 
Y ' ar^ 

2-a 


a 


2(1-a) 2(1-a) 


1 + 


4Kq‘^(l — a) 


(17) 


(18) 


Using, the Eq. CHI), we compute the integral (f Vr^drJ (see the Appendix) and use it in the Eq. (fTHll to find the 
expression of 


faJoo > a > 1 case: 
Using Eq. (IA.3I) . we get, 


g2^(r) ^ ^ 


ar 


6K(a — 1 ) 
2M 
f 

r4>\ 


^ — 1 ) ^ 


Q,i/4(4q2)3/4 / ^ f (AK,q^(a-l)) 

-F arcsm — 


1/4 ^ 


3/v1U(q: — l)i/4r I 




-1 


(19) 


where, F((j)\m) = /^^[l — TOsin^ 0] is the incomplete elliptic integral of the first kind. Thus we find all the metric 

functions. Here, the radius-squared of the 2-sphere at each value of the radial coordinate f is given by r^(f) = U(f)f^. 
We note that, for 1 < a < 2, at f = (areg^)^/^, the radius-square (r^) becomes zero. But, for a > 2, there is a non-zero 
minimum value of r^, i.e. r^(fo) = ^ 2 (a^i) ’ '"^here vq = [4Kg^(l — l/a)\^^^. This means, for 1 < a < 2, the spacetime 
is due to a point charge (q) of total mass M, but, for a > 2, it is distributed over a spherical shell in space. One can 
verify that, for a very small value of k, the solution converges to the known Reissner-Nordstrdm solution. The reason 
for this is that, as k becomes very small with respect to the value of a, becomes large, which then is the Maxwell 
limit for the electric field and the GR limit in the gravity sector. Eurther, we have examined the singularities of the 
spacetime and found that curvature scalar (TV) diverges at the position of the charge. Moreover, the existence of the 
horizon depends upon the values of the parameters a, k, q, and M. Some of these features are demonstrated in Fig. [T] 
and Fig. [H where the plots of the metric functions in the Schwarzschild coordinates are shown. 

In the Schwarzschild gauge, the line element [eq. ([6])] becomes 


ds^ = —gudt^ + grrdr^ 
where, gtt = t/(r)e^’^^’'^ 


(de^ 

U ( f )^_2 


• Sin' 


9d(t)^) 


grr = -rri-±e and, 


U(r) 


^ a a 
^-2 + 2 


i-b 




1/2 


( 20 ) 


(iv)—oo < a <1 case : 

In this case, using Eq. (|A.7|) we get. 


^2ip{r) _ 


= 1 - 


2M 


ar 


6k( 1 — a) 


^ ^ 4Kq^(l-a) _ ^ 


4g^ 


115 4Kg^(l — a) 
4’2’4’ ^^74 


( 21 ) 
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where, 2 -Fi(a, &; c; z) is the hypergeometric function with usual meaning. Here, for 0 < a < 1, the r^(ro) = 0, 
for fg = But, for a < 0 {i.e. for k < 0), there is a minimum value of the radius-square, which is 

r‘^{fo) = ( ] , for rg = 0. Therefore, for a < 0 or k < 0, the charge (g) and the mass (M) are distributed over 

a 2-sphere, but, for 1 > a > 0, it is a point charge. This solution also reduces to the Reissner-Nordstrom solution for 
very small «(«—>■ ±0). Here also, the spacetime is singular at the location of the charge and it may be a black hole 
or naked singularity depending upon the parameter values. The metric functions for this class of solutions are plotted 
in the Schwarzschild coordinates [eq. (EUl) ] for some parameter values and these are shown in Fig. [2] and Fig. [31 


B. Deviation from the solution with Maxwell field 


For K > 0, the solution with the Maxwell field (a —oo) is given by (using Eqs. |17I18I19I3D] '). 




-dr^ 


where, hi{r) = 1 - — - 


{r‘^ — Kq'^) {r‘^ + Kq'^)hi{r) 

2Mr 


-b 


{d0^ 

( 




F arcsin 


sm* 9d(j}^) 


^/r'^ + nq^ 


-1 


( 22 ) 

(23) 


Since, the charge is distributed over a 2-sphere of radius rg = (reg^)^/^, therefore r > {nq^Y/'^. For k < 0, the solution 
with the Maxwell field (a —>■ —oo) is given by the Eq. (I33|) . but with the replacement of hi{r) by h 2 {r), where 


^ 2(0 = 1 - x -2 - 


2Mr 


2^2 


-b 


Aq^r 


115 dKg^r^ 

2-^1 I 


3r2 _|_ ^q2 3(r4-b Kg^) \4’2’4’(r"^-b Kg^)^ 


(24) 


The deviation of the EiBI solutions with Born-Infeld electric field from that with the Maxwell field are demonstrated 
in Eig.|T]and Fig.|2]for k positive and negative respectively. In each figure, k- value is fixed, but a = takes different 
values for different solid lines. Therefore, Born-Infeld electric field parameter (6^) is varied without disturbing the 
gravity sector. All plots show the variation of metric functions in the Schwarzschild coordinates with physical radial 
distance r. The plots are obtained using Eqs. [T71 [HI 11912111331130113311331133] . For a given value of a, an appropriate 
aforementioned equation has been used. In some plots r starts with a nonzero minimum value which is due to the 
charge s dist ributed over a 2-sphere of finite radius. 

In Fig. [l(a)| we clearly see that the metric function [gu) for different values of a = 4 k 5^ deviate from the case with a 
Maxwell field (black dashed line) at small r. However, at large r, plots for different a approach the black dashed line 
(Maxwell) coinciding with it asymptotically. We also see that, for larger values of a or 6^, the amount of deviation 
reduces. In Fig. 1(a) and Fig. [1(b)] we note that for the chosen parameter values {k = 1.0, g = 0.5, M = 1.0) the 
deviatio n in t he horizon radius from the case with a Maxwell field is negligible. We note a clear deviation ( Fig. 1(c) 
and Fig. [l(d)| ) in the horizon radius for a different value of g = 0.9. Here, the horizon radius increases with decreasing 
a, or 6^. In Fig. [T] one can see that for a = 0.5 and a = 1.5 the metric functions (both gtt and grr) begin from r = 0 
indicating a point charge. In contrast, for a = 3.0 and the Maxwell limit (black dashed line), the metric functions 
are limited to a minimum radius of the 2-sphere where the charge is distributed. From the Fig. [31 one can also note 
the existence of a minimum radial distance for all a with k negative. 


C. Deviation from the Geonic blackhole solution 


Additionally, we note that, in the limits- a —?► ±0 and k —>■ ±0, but for a finite 6^ we get back the geonic black hole 
solution. The geonic blackhole is a solution in GR coupled to a Born-Infeld electric field due to a point charge [43 |. 
In a geonic black hole scenario, a distant observer associates a total mass which comprises M (the black hole mass) 
and a pure electromagnetic mass stored as the self energy in the electromagnetic field. If M is zero, the spacetime 
becomes regular everywhere. 

The line element for a geonic blackhole is given as (from Eg. 1311) . 


ds^ = 


-ge{r)df 

2M 


dr^ 

9e{r) 


\de^ 


2g2 


sin^ OdcY) 
4g^ 

-1 —— ' 

3r2 ■ 


Fi 


1 1 5 


(25) 


where, geY) = 1 


r 


3 




(26) 
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r 


(a)«; = 1.0,1? = 0.5, M = 1.0 


(b)K = 1.0, (? = 0.5, M = 1.0 




r 

(c)k = 1.0,1? = 0.9, M = 1.0 


r 

(d)K = 1.0, g = 0.9, M = 1.0 


FIG. 1. Metric functions in Schwarzschild coordinates are plotted and compared with the EiBI solution for the Maxwell electric 
field (black dashed line), gtt and Qrr are the functions given in the equation (I20II . In the plots, k = 1.0, while a = takes 
different values, which means the solutions are with different values of control parameter (b^) for the nonlinear Born-Infeld 
electric field. 


where, r is again the physical radial distance. Here, r can take any positive value from zero to infinity. The deviation 
of the EiBI solutions (with different k) from the geonic blackhole solution (the GR limit) is demonstrated in Fig. [3] 
with a fixed value of the Born-Infeld electric field parameter (6^). Here, the gravity sector is varied by changing the 
value of K without disturbing the matter sector. 


The deviation in the metric functions from the geonic blackhole solution (brown dashed line, k = 0) is prominent for 
different values of k in Fig. However, at larger r, all plots approach the brown dashed line. For a larger magnitude 
of K, the deviation is more. In Fig. |3(a)| and Fig. |3(b)| we see that for the given parameter values, the double horizon 
exists in the case of the geonic blackhole, though it is absent for other solutions with different values of k. From 
Fig. 3(b) 3(d)[ and |3(f)[ we note a generic feature that for k positive the horizon radius decreases with increase in 
magnitude of k. But, it is the opposite for negative k. 


IV. THE n = 1 CASE 

A. The line element 


For 0 = 1, the quadratic equation [Eq. (flO)) ] becomes linear in V. Hence, the metric functions V, U and 
[Eqs. (fT71) . (fT3]) . (IT^ ] become 



U = l 


Kq 


and 


= 2 ^ = 1 _ 


2M 


(27) 
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(a)K = -1.0, q = 0.5, M = 1.0 


(b)K = -1.0, (j = 0.5, M = 1.0 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 


r 

(c)k = -1.0, g = 0.9, M = 1.0 


r 

(d)K = -1.0, g = 0.9, M = 1.0 


FIG. 2. Metric functions in Schwarzschild coordinates are plotted and compared with the EiBI solution for the Maxwell electric 
field (black dashed line), gtt and grr are the functions given in the equation Here, in these plots also, K-value is fixed but 
negative and consequently a is also negative. 


Then, the physical line element becomes, 


= — 






{de^ + sin^ 0d(j)^) 


(28) 


Here, the radius-square: r^(r) = becomes zero at f = (kq^)^/^. So, this describes a point charge q with the 

total mass M. This line element [Eq. (E51) ] clearly shows that, for k = 0, it becomes the Reissner-Nordstrom solution. 
The auxiliary line element is also the same as Reissner-Nordstrom. 

The source for such a spacetime has the following form of the stress-energy tensor in the framework of EiBI gravity: 


= diag. \ - 


8Tr{r^ — nq^) ’ 87r(r^ — nq^) ’ + nq^) ’ 87r(r^ -|- nq^) j 


One may verify that the stress energy satisfies the weak/null energy conditions [^ . 
In the Schwarzschild gauge, the line element takes the form 


ds^ = —f{r)dt^ -1- 


2H 


{r‘^ + \/r^ + f{r) 


dr'^ + r^ (d6»2 -bsin^ Odcj)^) 


(29) 


(30) 


where, the redshift function, /(r) 


( 1 1 w \ 

-| 2V2M 1 2q^ 

1 J 

^r'^+y/r'^+Anq'^ r‘^+^/r^+'^Kq^ 


At the horizon- 
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(a)fe2 = 1.0,1} = 0.9, M = 1.0 



r 

(b)fe2 = 1.0, g = 0.9, Af = 1.0 


- Geonic Blackhole 

- K= 1.0 

_ k = 2.0 

- /f=5.0 

K= 10.0 

1 




15 .. 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 


r 

(c)b2 = 1.0, g = 0.5, M = 1.0 


r 

{d)b^ = 1.0, g = 0.5, Af = 1.0 




r 

(e)b2 = 1.0, g = 0.7, M = 1.0 


r 

(f)b2 = 1.0, g = 0.7,= 1.0 


FIG. 3. Metric functions in Schwarzschild coordinates are plotted and compared with the geonic blackhole solution (brown 
dashed line), gu and Qrr are the functions given in the equation (I20l). Here, in these plots, the Born-Infeld electric held 
parameter (&^) has the hxed value = 1, but EiBI gravity parameter k is varied. 


radius, rhz, fi^hz) = 0. The expression for the horizon-radius is therefore given by. 


' hz 


(^M -q^y - 

(M±y/M^ - q^y 


Kq 


(31) 


Here, we note that a horizon does not exist for q > M, similar to the Reissner-Nordstrom case. For q = M, we have 
= q'^ — n- So, unlike the Reissner-Nordstrom extremal solution, for each value of g = M, there is an upper limit 
on K for the existence of the horizon, (see the left panel of the Fig. 0]). From the right panel of the Fig. 01 we note 
that, for q < M, for sufficiently small values of k there may be double horizon. For a range of higher values of k a 
single-horizon exists, and for even more high values of k the horizon vanishes. However, in the Reissner-Nordstrom 
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counterpart, we always have the double horizon. 


1.0 


0.8 


0.6 - 


0.4 - 


0.2 


0.0 


=\l A K ,q = 'M. 


Black4iole 


0.0 


0.2 


0.4 0.6 

q 

(a) 


0.8 


3.0 
2.5 
2.0 
^ 1.5 
1.0 
0.5 


1.0 


0 . 0 ^ 

0.70 


1)2->1/(4k) 


Single horizon 




\ O - 




17 . 0 ' 


0.75 0.80 0.85 0.90 

qlM 
(b) 


0.95 1.00 


FIG. 4. (Left panel)the parameter space: — k showing the existence/non-existence of horizon, for q = M. (Right panel) the 

same plot for q ^ M case. 


One can verify that, in general, all the scalar quantities, like the Ricci scalar (TZ = g'^^Rij), and the 

Kretschmann scalar diverge at the location of the point charge, i.e. at f = or, r = 0. However, 

for a special choice, q = i/k/ 3 and M = 2i/k/ 3-\/3, the Ricci scalar becomes regular everywhere. But, the spacetime 
still remains singular at the location of the point charge, for such a choice, since the Kretschmann scalar and RijR^^ 
diverge. Fig. [5] demonstrates the singularity for two different value of the parameters. So, the spacetime is eventually 
singular, but may or may not be covered by a horizon. However, we show that the strength of the electric field 
at the location of the charge remains finite. For this, we consider a static observer having the four velocity, = 
0, 0,0}. The electric field vector is defined as So, the square of the strength of the electric 

field is E'^ = E^E^ = At the location of the charge, E"^ = /b 2 = 6^ = 1/4k (in this case, 


f{ro) = 0). The energy density observed by such a static observer is pobs = 


i,2 f ^y2p4+g2/b2 ^ 

4ti- \ FP J. 


In 


this case, Pots = STT(f^-K,q^) ■ diverges at the location of the point charge (r = tq = We note that 

the non-singular feature of the strength of the electric field remains even in the solution of EiBI gravity coupled to 
Maxwell’s electrodynamics [1, and this is due to the Born-Infeld like structure in the gravity sector. However, 
this does not ensure the nonsingularity of the spacetime. In im, the charge is not point-like, but distributed over 
a spherical surface. Pobs remains finite, though the scalar invariants diverge at this surface. But, in our solution, 
both the energy density and the curvature scalar are singular. A non-linear correlation between the gravity and the 
physical matter sector in the EiBI theory is thus revealed. 


B. Null geodesics 


Let us now turn to the null geodesics in this spacetime. On the equatorial plane {6 = 7r/2), for null geodesics, we 
have, from Eq. 


^ 1 - 


Kq 


-b 


2r^(r) 


1-b 


Kq 


1 - 


2M 


q 


'Y 


(32) 


where, E = {1 


d(p • _ ^ 

dA ’ ' ~ df dX^ 


^ i, L = are two conserved quantities and i, (p, r are defined a,si = p 


A is the affine parameter. So, from Eq. (l32)l we identify the effective potential for the null geodesics as 


, , / Kq^\ f 2M 

= + (i- — 


(33) 


where, = (r^ + y'Y+AKcY)I■ The aforementioned effective potential is useful to study the propagation of 
gravitational signals in the current scenario. We discuss photon propagation in the following subsection. 
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q =M = 0.4,/(: = 1.0 


q = 0,9,M = 1.0,/f = 0.5 




(a) 


(b) 


FIG. 5. (Left panel)the ricci scalar {TZ) is plotted as a function of r for the parameter values, q = M — 0.4, k = 1.0. The plot 
shows that TZ diverges at the location of the charge, i.e. f = or, r = 0. (Right panel) the same plot for a. q ^ M case. 

Here, TZ diverges to negative values. 


C. The propagation of photons 


Photons (lightlike particles) in Born-Infeld electromagnetism (BI photon), no longer propagate along the null 
geodesics of the background spacetime geometry. This is generally true for nonlinear electrodynamics [s^. Instead, 
the BI photon trajectories are determined by the null geodesics of the so-called effective geometry, given by 


9eff 




(34) 


where, F = |Fq,^F“^ and defines the metric of the background spacetime. Such an effective geometry is induced 
due to nonlinear self-interaction of the photon. Here, in Eq. (1341) . we assume the same F^j^ which generates background 
spacetime geometry. Then, the effective line element becomes 


dslff = - ( 1 


nq 


1 - 


2M 


dt^ 


(i + ^) (i + ^)^ 


- —^dr + -7 -<- 

+ ( 1 -^) 


{dO'^ + sin^ 0d<jy^^ (35) 


Since our interest is to study the null geodesics of the effective geometry, we take out the conformal factor 
and rewrite the effective line element as. 


dslff - - ( 1 - 


2M 

f 


dt^ 


( 1 -^ + $)'^ ( 1 -^) 


(d0^ 


Sin' 


' 0di^^) 


(36) 


The effective geometry becomes exactly the same as the Reissner-Nordstrdm solution when k = 0, which is expected. 

, shows that the area radius, 


= r' 


f f‘^ + K.q^ \ 


The analysis of the squared area radius in the effective geometry, 

or, the physical radial distance r becomes infinity both at r = Fq = (zcg^)^/^ and f —>■ 00 , but has a minimum 

~ ~ ( -ji 2 ) (similar to the throat radius square of a wormhole). This implies that 

for each physical radial distance r G [rth,oo], there exists a r G [fth, 00 ] and a f G [ro,fth]- But these two sets of ranges 
of r describe two different geometries. We are interested in studying the propagation of photons coming from and 
going towards asymptotically flat regions of the given spacetime. Hence, we focus on rth < f < 00 , which describes 
an asymptotically flat geometry. 

Now, for the null geodesics, we have 


1 -2 
—r 
2 


2f2 


1 - 


2M 


r 

jt2 


-4 2 

r — nq 


F2 


(37) 
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where, E = 1 - ^ 


become 


r)t 


and L = r^{r)(j) are two conserved quantities. One can also rewrite the Eq. (1371) in terms 

2(^\ — ^2 f f‘^+Kq^ ^ 




1 - 


f2 ( 


r 


2M 

q^' 


+ 


So, the effective potential for the BI photons 

(38) 


-4 2 

r — Kq 




i^/M^ = 0.0,q/M = l 



k! M ^ = 0.9,q/M = l 



k/M ^ = 6.0,q/M = l 



r/M 


r/M 


r/M 



(d) 


(e) 


(f) 


FIG. 6. (Top panel)Tlie effective potential [Eq. (I38II ] is plotted as a function of the dimensionless quantity r/M, where r(f) is 
the radius of the 2-sphere labelled by the radial coordinate f and M is the mass. In the parametric plot, the (f) values are in 
the range rth < r < oo. In the plot, Veff = 2M^I4///I/^. The parameter values are specified on the top of the figures. The 
vertical red dashed line in each plot shows the correspondence between the photon sphere radius and the local maxima of the 
effective potential. (Bottom panel)The exact deflection angle {Arj)) is plotted as a function of rtp/M, where rtp is the physical 
radial distance at the turning point. In each figure, A(j) rises sharply along the vertical red dashed line. The corresponding 
rtp/M relates to the photon sphere radius. 


We plot the effective potential for different values of k/M'^ (a dimensionless parameter) and extract the essential 
features. A few plots are shown in the top panel of the Fig. [SI In these plots, the effective potential is rescaled as 
Veff = 2M^Veff /L? and is plotted as a function of the dimensionless quantity r/M. In these parametric plots, the 
parameter (r) value ranges from fth to infinity. The maxima of the effective potential correspond to the unstable 
circular orbits for photons. The corresponding spherical surface is called the ‘photon sphere’ [13, [13 [11|. In the 
Reissner-Nordstrom solution, the radius of the photon sphere depends on q and M. In our solution, we note that, it 
depends additionally on k. As an example, we consider the case for q/M = 1. In the Reissner-Nordstrom solution, 
the radius of this photon sphere is rpg = 2M = 2q (see Fig. 6(a)). In our solution, we have a different picture. Here, 
from the figures of the top panel of Fig. |6l we see that as n is increased the radius of the photon sphere increases. 
This feature becomes more clear when we look at the plot of deflection angles {A(j)) in the bottom panel of Fig. [51 
Near photon sphere, the deflection angle becomes infinitely large, indicating multiple full rounds of light rays around 
the lensing object. Consequently, a sequence of a large number of highly demagnified relativistic images may arise 

Analysing the plots of the effective potential for greater values of k/M^, one can easily verify that the radius of the 
photon sphere {rps) increases but becomes closer and closer to the throat radius (rth) as ^/M^-value increases. 
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Strong deflection: In this effective geometry, the exact expression for the deflection angle 
A(j) = 2|(/i(oo) - </>(rtp)| - TT 


becomes, 


= 2 




f^(f^ + Kq^) 


' tp 


rfp - 

rfp + 


1 - 


2M 

rtp 


' tp 


1 


-4 2 

r — K.q 

+ Kq^ 


1 - 


2M 


- 1/2 


dr 


(39) 


=2 


where, r^{rtp) = . Here, rtp indicates the turning point, i.e. |^|rtp = ^\rtp = 0. The deflection angle 

is plotted as a function of the dimensionless parameter vtpfM and the flgures are shown in the bottom panel of 
Fig. m Fig. \6{d)\ (for k/M^ = 0, q/M = 1) shows that, in the extremal Reissner-Nordstrom case, the deflection 
angle increases rapidly (>> 27r) near the photon sphere {rtp = rps = 2M). From the next two figures (Fig. 6(e) and 
Fig.M, we see that the nature of plots do not change. However, as the k/M"^ value increases, it is clearly seen that 
radius of the photon sphere increases from the value 2M. 


Weak deflection: Let us now look for the leading order contribution of k in the weak-deflection angle formula. We 
compare this result with the Reissner-Nordstrom case. We obtain the weak deflection angle perturbatively M- For 
a spherically symmetric static line element, ds^ = —A{f)dt^ + B{f)df^ -I- f^C{f)dn^, the radial geodesic equation is 


2Bf + B'P + A'r - = 0 


(40) 


With the conserved quantities, E = At, L = Cr'^cp, and the new variable u = 4, for the null geodesics, the Eq. (HOI) 
becomes 


d'^u C I pd fC 
dcjfl B^ ~~2^ du \B 


1 d f C 


26g du V AB 


(41) 


where, 6o = LjE is the impact parameter. In the Reissner-Nordstrom case, (^A = B ^ = 1 — 2Mu + q^u^, (7 = 1), 
the Eq. (HD) becomes. 


d?u 

d(jfl 


+ u = 3Mu — 2q u' 


2„,3 


(42) 


In our case, A = B ^ = (l — 2Mu + q'^u'^), and C = {\ + Kq^u^)/{1 — nq^u'^). We assume small n {nq^u^ « 1) and 
write the Eq. dun keeping only the first order terms in k, 


d?"h 

d(jy 


+ u = 3Mu^ —[ 2 — 


8k 


q^u'^ - Gnq^u^ -|- MKq^Mu*^ - SKq^u 




(43) 


We set e = Mun, where un = l/?'Ar and bo = r^ = r^ ({cn is the Newtonian distance of closest approach). 

\'^N J 

Further, we define a dimensionless variable ^ = u/um- Thereafter, we write the Eq. (H51) to second order in e. 


0 + 


2 - 


8k 


^^3 


(44) 


Expressing ^ as ^ + e C 2 +., and substituting them in the Eq. (lUl) we find the following set of equations, 

by collecting terms of different order: 


d(jfl 

dcffl 

dcjfl 


+ Co =0 


+a - 3 a' = 0 


+ a ~ 6 aa + 2(1 — 


4h; 


a' = 0 


bl 1 


(45) 

(46) 

(47) 


Solving the equations (1^ . (l46)) . and (iTf)) we get the approximate inverse radial distance, 


u ~ Un coscj) + Mun { 2 ~ 2 cos(20) ) -|- 
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20-4(1-—)^ 

^ bl’M^ 


^sm(j)+ ( 1 -h ^1^(1 - ^) ) cos(30) 


(48) 
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Assuming, u —>■ 0 at (/> = 7r/2 + 5, we solve for i5 using the Eq. (I48|) . For a small deflection angle, we find 


,5 


2Mun + 


8 


5-(l 


bl'M^ 


0{M^u%) 


(49) 


The deflection angle is A(j> = 26. To get the first order in k contribution in the deflection angle, we use bg « I/mat. 
At the distance of closest approach, rtp = l/u(0). Converting from un to m( 0), substituting u{0) = 1/ftp, with 
ftp ~ ftp (l ~ /ftp + , we finally obtain the approximate weak deflection angle formula as. 


A(j)^ 


eak 


AM 

‘^tp 



(50) 


In this formula, the term involving k is the new contribution to weak deflection angle. If k = 0, the expression becomes 
the same as in the Reissner-Nordstrom case [^ . By extending the calculation upto the fourth order perturbation 
(i.e. one can verify that the sign of the leading order in k term that appears in the fourth order correction 

(1/r^p-term) in the deflection angle, is still positive. This term involving k increases the net deflection angle. As the 
value of K is increased, the deflection angle also increases further. Earlier, we have seen that, with the increase of k, 
the radius of photon-sphere increases. Thus, the study of both strong and weak deflection seem to suggest an overall 
attractive effect of k. 


V. EFFECTIVE POTENTIAL AND PHOTON PROPAGATION FOR a ^ 1 SOLUTIONS 


In this section we investigate the effective potential, photon sphere, and light propagation for a ^ 1 solutions 
through numerical graph and plots. We also compare these with that for the Maxwell electric field by taking different 
values of (or, a with fixed k). 

The effective metric for the BI photon propagating in the background spacetime which is the EiBI solution with 
Born-Infeld elecric field (a 1), is derived from Eq. (IMl) using Eqs. [51 [7] as, 

dslff = —U{f)e‘^'^^^'>dt‘^ + U{f)e~‘^'^^^'>df‘^+ (^ {d6^ + sin^ 6d(j)^) (51) 

jj y aV (r)r^ J 


where, t/(r), Vir), and have the expressions given by the Eqs. [171 HU and [19] or, [21]. Using Eq. (EID, and 

following the same analysis as described in section IV, the effective potential for BI photon (for a 1) is derived as. 


Veffif) 


where, r^(r) 


2r^(f) 

/ aV‘^{f)f'^ + Anq^ \ 

\ aV{f)f‘^ ) ' 


(52) 

(53) 


Here, the effective potential as a function of the physical radial distance, i.e. Veff{f), is to be determined in a 
parametric form as given by Eqs. [([SI]), dSSl)]. 

On the other hand, from Eq. (I22L we get the effective potential for the Maxwell photon as 


eff 


(r) 


r^hi(r)L^ 

2(M — nq^) 


(54) 


where, hi(r) is given by Eq. (ITS)) . In Eq. (1571) . the effective potential is written directly as a function of the physical 
radial distance r. For k < 0, hi{r) is replaced by h 2 {r) (given by Eq. (l24l) l. In Fig. [71 variation of effective potential 
for BI photons as compared to Maxwellian photon are shown. For plotting, we use Eqs. |54I521 ISS) . Similar to the 
a = 1 special case, we consider f € [fth, oo], where, = yjKcp-/a/l + a{'/2— 1)], for a range of a. From the plots in 
Fig. [71 it is evident that the radius of the photon sphere {fps) for BI photons is smaller than that of Maxwell photons 
(irrespective of positive or, negative k). Furthermore, the photon sphere decreases as the values of b^ decreases. 

Let us now write the expression for the deflection of light (BI photon) from the null geodesic equation for the 
effective metric given by Eq. (EH). We get the expression for the deflection angle as. 


A(j) = 


2 


/■- U{f) 



f‘^{ftp) f'^{f) 


- 1/2 


dr — TT 


(55) 
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0.04 

0.02 

^ 0.00 
- 0.02 
-0.04 


FIG. 7. Effective potential (14//) for BI photons are plotted and compared with that of the Maxwellian photon (black dashed 
line) propagating in the background of EiBI spacetimes. In both plots, k value is same (but of opposite sign). However, 
a — AKb^ takes different vaues in the plots. 



where, r^(r) is the physical radial distance square given by Eq. (15^ and rtp is the coordinate value corresponding to 
the physical radial distance at the turning point (rtp). Again, for the Maxwell photon, we get the deflection angle. 


= 2 


fOO ^ 

rfphi(rtp) 

r'^hi(r) 

frtp + Kq^ 

1 

1 

to 

4 2 


-I -1/2 


dr — ■ 


(56) 


where (1561) is expressed directly in terms of physical radial distance (r). In Fig. [H the variation of the deflection 
angle of light as function of rtp is shown for different a values and it is compared with that for the Maxwell field. 
To visualize the correspondence between the effective potential for photons with the deflection angle of light, same 
parameter values are used in the plots. For plotting, Eqs. [551 ESI ES] used. 




Hp Hp 

(a)K = 1.0, L = 1.0,1} = 0.7, M = 1.0 (b)K = -1.0, L = 1.0, g = 0.9, M = 1.0 

FIG. 8. Deflection angle (At/i) for BI photons are plotted as a function of the physical radial distance at the turning point (rtp) 
and compared with that of the Maxwellian photon (black dashed line) propagating in the background of EiBI spacetimes. The 
deflection angle diverges rapidly when the turning-point-radius (rtp) becomes closer to the radius of the photon sphere (rpa). 
In both plots, the parameter values are the same as in Fig. (3 

From the plots of both the effective potential (Fig. | 3 ) and strong deflection angle (Fig. [ 5 ]), we note that with the 
increase of the value of 6^, the radius of the photon sphere (rps) decreases. 
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VI. CONCLUSIONS 


In this article, we have derived new classes of spherically symmetric static solutions in EiBI gravity coupled to 
Born-Infeld electrodynamics. The spacetimes demonstrate the effect of coupling a non-linear electrodynamics (Born- 
Infeld electrodynamics) to a modihed theory of gravity (the EiBI gravity). All solutions are analytical. We summarize 
our results below, point wise: 

(i) We have used a non-standard choice for the physical and auxiliary metrics to find the spacetimes in EiBI gravity. 
The standard ansatze for the physical and auxiliary line elements as employed in [1, Idcil - I^ , led to a complicated 
set of coupled differential equations. We note that our non-standard choice is far more useful insofar as solution 
construction is concerned. Similar unconventional choices may be used in other cases with the hope of simplifying 
the set of equations we need to solve. 

(a) We obtain two classes of spacetimes in section III. All the known results for different electrically charged, 
spherically symmetric, static scenarios can be obtained by taking different limits of a = and k. In general, the 
spacetimes we obtain are singular at the location of the charge, where the charge can be localized about a point or 
distributed over a spherical shell (2-sphere). The singularity may or may not be covered by an event horizon, i.e. it 
may be a black hole or a naked singularity. However, all solutions converge to the Reissner-Nordstrom solution at 
large distances. 

{in) The effects of modifications in both gravity and matter sectors are demonstrated through the analysis of the 
metric functions for fixed 6^, varying k and vice-versa. In our analysis, we find some interesting features. For a fixed 
K, the horizon radius increases with the decrease of 5^ for some parameter values {K,q,M). If 6^ is fixed, the horizon 
radius decreases with increasing k {k > 0). 

(iv) Although the results are analytical, the expressions for the metric functions are complicated. However, a special 
choice, a = 1 (k > 0), leads to a simple solution that looks similar to the Reissner-Nordstrom. We analyze this special 
solution in detail in section IV. We find that this spacetime also possesses the same generic features mentioned above 
in (a). Subsequently, for the a = 1 line element, we investigate the null geodesics in the effective geometry for BI 
photons, as well as gravitational lensing. We find that as k is increased, the radius of the photon sphere increases. 
We also find the leading order contribution of k in the expression for the weak deflection angle. The overall sign of 
this term is positive but it depends on the value of k. 

(u)In section V we investigate the effective potential, photon sphere, and strong deflection angle for a yf 1, through 
numerics and plots. We compare with results for a Maxwellian electric field. Analyzing the plots, we find that for a 
fix value of k (irrespective of its sign), the radius of the photon sphere decreases as the value of 6^ is increased. 

In summary, using non-standard ansatze for the physical and auxiliary metrics, we have been able to obtain analytical 
solutions in Born-Infeld gravity coupled to Born-Infeld electrodynamics. One may explore strong gravitational lensing 
further using the approach in (^ . It is possible that there are other non-singular solutions as well which are non¬ 
singular from both the gravitational and electrodynamic perspectives. For this, it is essential to scan the (k, 6^, ■^) 
parameter space more thoroughly. In , the author has shown that in the case of EiBI solutions with a Maxwellian 
electric held, there exists a region in the parameter space (k, q, M) which allows the construction of a Lorentzian 
wormhole without violating the weak or null energy conditions. In our solution, a similar analysis is possible. One 
can also study the black hole thermodyanmics using the solutions presented here. We hope to work on these issues 
and communicate our results in the near future. 
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Appendix: Integral results 

Here, we compute the integrals / V(f)f^df, in detail, which have been used in the section HIIl to obtain the analytical 
expressions of two different classes of spacetimes [Eq. (IT^ and Eq. (I?T]) ] . These integrals are non-trivial and may not 
be computed correctly otherwise (e.g. using Mathematica). 


17 


(a)oo > a > 1 case: 


Vr^dr = — 


(2 — a)r' 


^3 


_ W(o-l)^- 


J 6 (a — 1) 2{a — 1) J \ a 

The last term on the R.H.S. of the Eg. lA.ll is a non-trivial integral which we evaluate below: 

2r'^df 


(A.l) 


— fddr = r-^/r^ — /3 — 


- P 

= f\/— jd ~ 2 y — fddr — 2(3 
2 / 33/4 


( Using partial integration) 
df 


- P 


j0=i' 


Now, 


dy 


du 


- 1 


V 1 — u4 ’ 


\/ 2 /'‘ - 1 
(u = l/jz) 


= -F (arcsin(l/y)|-l) + C^xt 

So, J \/f4 - (3dr = if\/f4 - /3 + ^arcsin(/3^/^/r) -ij - Cext) 


(A.2) 


where, F{(j>\m) = /q^[1 — msir? 9]~^^‘^d6 is the incomplete elliptic integral of the first kind. Here, (3 = 

(/3 > 0) and Cext is an additional constant term appearing due to the fact that Elliptic integrals are defined with the 
zero lower limit and hence, these are not the so-called indefinite integrals. The Eq. (lA.ll) becomes. 




Vr dr = — -I- 


3 6 (a-l) 

a4/4 (4^92)3/4 

3(a- 1)1/4 


/ _ - 1) _ 1 
V ar^ 

{AKq^{a — 1))^^^ 


F ^arcsin 


q,1/4j. 


-1 - C, . 


(A.3) 


Now, to fix up the extra constant C^xt^ we use the fact that, in the limit of a —> 00 and k —> 0, the full expression of 
must take form of the Reissner-Nordstrom solution (i.e. e^’/’ —>• 1 — Under this constraint, Cgxt = 0. 


(b)—oo < a < 1 


Vr^dr = — 


(2 — a)r^ 


■ J \/r^ + 'ydr 


6(a- 1) 2(a- 1) 

where, 7 = 4«9 (7 > 0). We evaluate the last integral on R.H.S. of the Eq. (jA.4ll below: 

J idf = irVf4 -h 7 -k J , y (r = 7^'^'‘2/) 


Now, 

So, 


dy 


, _= = —\/iF {i arcsinh(Viy) — 1 ) -b C'f 

y/ipi V \ yj ) 

^ViF (i arcsinh{Vir^~^^^) — 1 ^ — 


\Jr^ F jdr = -r\/ r4 -|- 7 — ^ 


where, i = So, the Eq. (IA.4I) becomes 


Vr dr = — 


r ar 


3 6 (q;- 1) 

ai/^(4Kg^)3/4 

3(l-a)i/4 


^ 4«;g3(l — a) 

1 H-=4- 1 

ar^ 


ViF yarcsinh 




{dKq^{l — a))i/4 


-1 - C'^xt 


(A.4) 


(A.5) 


(A. 6 ) 
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Similar to the earlier case, we fix up the extra constant by taking the limit of a —>• —oo and k —> —0. 

Then also, the full expression of takes the form of the Reissner-Nordstrom solution. Under this constraint, 


CLt = -fr 


dy 


\/y‘‘ + i 


= — ■ Eq. (IA.6I) can be rewritten further as, 


Vr^dr = — 


3 6(q;-1) 


W(l - a) ^ 
af'^ 


f 4Kq" 


V 3r 


115 Anq'^il-a) 
^^'4’2’4’ af4 


(A.7) 


where, we use 


ViF ^farcsinh (^Viy^ —1^ 


r^(i/4) 

4y/Tr 


dz 


\/i4TT 


Here, 2 ^ 1 ( 0 , b; c; z) is the hypergeometric function with usual meaning. 
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